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THE MID-INFRARED EMISSION OF LOCAL LUMINOUS INFRARED 

GALAXIES 

Almudena Alonso-Herrero 1 



O 
O 

(N 

O 

Q 

o 



> 

in 

(N 
O 

a: 

6 

in 
c3 



RESUMEN 

En este artfculo se combinan la sensibilidad de Spitzer y la alta resolution espacial de T-ReCS en el Gemini- 
South para estudiar las propiedades en el infrarrojo (IR) medio de galaxias luminosas infrarrojas (LIRGs) 
locales (d < 75Mpc). Las imagenes de T-ReCS en 8 — 10 ^m nos han permitido separar la emision nuclear 
(formacion estelar y/o AGN) de la de las regiones Hn en las zonas centrales (3 — 7kpc). De la comparacion 
de los cocientes 8/xm/Paa integrados y de regiones Hn resueltas se deduce la existencia de una componente 
de emision difusa en 8 /im. Esta emision difusa no esta directamente relacionada con las poblaciones estelares 
ionizantes, y puede ser tan luminosa como la que proviene de regiones Hn. Por lo tanto, aunque la emision en 
el IR medio es un buen indicador de la tasa de formacion estelar (SFR) de galaxias con polvo, las calibracioncs 
han de hacerse usando las propiedades integradas dc galaxias, y no las de regiones Hn unicamente. Se presenta 
una calibration dc la SFR en funcion de la luminosidad en 24 fim que es aplicable a galaxias lejanas con polvo. 

ABSTRACT 

In this paper we use the complementary imaging capabilities of Spitzer (sensitivity) and Gemini-South/T- 
RcCS (spatial resolution) to study the mid-IR properties of local (d < 75 Mpc) LIRGs. The T-ReCS 8 — 10 (im 
imaging observations of LIRGs have allowed us to spatially resolve the nuclear emission (star formation and/or 
AGN) and that of H n regions in the central 3 — 7 kpc regions of LIRGs. From the comparison of the 8 ^m/Paa 
ratios of the integrated vs. resolved H II regions of LIRGs, wc infer the existence of an 8 /im diffuse component, 
not directly related to the ionizing stars, that can be as luminous as that from the resolved Hn regions. We 
conclude that although the mid-IR integrated luminosity of galaxies undergoing dusty, intense star formation is 
a good indicator of the star formation rate (SFR), the empirical calibrations should be based on the integrated 
emission of nearby galaxies, not that of H n regions alone. To this end we provide a calibration of the SFR in 
terms of the integrated 24 ^m luminosity that can be used for distant dusty galaxies. 
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1. INTRODUCTION 

The importance of infrared (IR) bright galaxies 
has been increasingly appreciated since their discov- 
ery more than 30 years ago (Rieke & Low 1972), 
and the detection of large numbers by IRAS. Al- 
though the Ultraluminous IR Galaxies (ULIRGs, 
Lm, > 10 12 L Q ) get much of the attention because 
they are so dramatic, Luminous IR Galaxies (LIRGs 
— -Lir = 10 11 to 10 12 L Q ) are much more common, 
accounting for ~ 5% of the local IR background com- 
pared with < 1% for the ULIRGs (Lagache et al. 
2005). LIRGs take part in the controversy over the 
formation of AGN, and its relation to star forma- 
tion and high levels of IR emission (see review by 
Sanders & Mirabel 1996), and should also include 
former ULIRGs where the SF is dying off. Deep 
Spitzer detections at 24 fim are dominated by LIRGs 
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(Le Floc'h et al. 2005; Perez-Gonzalez et al. 2005), 
and LIRGs contribute nearly 50% of the cosmic IR 
background at z ~ 1 (Lagache et al. 2005). Interpre- 
tation of the Spitzer observations at high-z therefore 
depends on a deep understanding of LIRGs in the 
local universe. 

With the advent of a new generation of mid- 
IR instruments with superior sensitivity and spa- 
tial resolution we can now investigate whether the 
mid-IR emission is a good indicator of the star for- 
mation rate (SFR) in dusty galaxies. The mid- 
IR emission has the advantage of not being af- 
fected by the contribution from cold dust heated by 
old stars that may dominate the far-IR luminosi- 
ties. Spitzer observations of nearby normal galax- 
ies show similar 8 — 24 (jxn and hydrogen recombi- 
nation line morphologies (e.g., Helou et al. 2004; 
Hinz ct al. 2004; Calzetti et al. 2005, CAL05 here- 
after; Perez-Gonzalez et al. 2006, PGPG06 hcrc- 
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after; Alonso-Herrero et al. 2006b). At longer wave- 
lengths (160 /im) the individual Hn regions (current 
star formation) contribute much less than at 24 and 
70 /zm (Hinz et al. 2004). From a quantitative point 
of view there is a good correlation between the Paa 
or Ha luminosity (corrected for extinction) and the 
24 /im luminosity of resolved Hn regions (CAL05; 
PGPG06). However there are galaxy-to-galaxy vari- 
ations in the 24 /im luminosity to SFR ratios of UV- 
selected starbursts and ULIRGs (CAL05). CAL05 
and PGPG06, on the other hand, questioned the use 
of the IRAC 8 /im emission as a SFR tracer, whereas 
Wu et al. (2005) find a good correlation between 
the 8 /zm (and 24 /zm) and Ha emission of local star 
forming galaxies in the Spitzer First Look Survey. 

We present a detailed study of the mid-IR prop- 
erties of local LIRGs taking advantage of both the 
sensitivity of Spitzer imaging observations and the 
high spatial resolution of the Thermal-Region Cam- 
era Spectrograph (T-ReCS; Telesco et al. 1998) on 
the 8.1m Gemini-South Telescope. In particular we 
investigate the issue of whether the mid-IR emission 
is a good tracer of the SFR in dusty galaxies. 

2. THE SAMPLE OF LOCAL LIRGS 

The sample of 30 LIRGs (34 galaxies) was cho- 
sen to be volume-limited to a velocity range of v — 
2750 - 5200 km s" 1 (d = 35 - 75Mpc, H = 75 km 
s" 1 Mpc^ 1 ) so that the Paa emission line (A rest = 
1.876 /im) could be observed with the fiST/NICMOS 
F190N narrow-band filter. We also obtained NIC- 
MOS continuum images at 1.1, 1.6, and 1.87/zm . 
For the velocity range used here our sample contains 
- 80% of all LIRGs in the IRAS Revised Bright 
Galaxy Sample (Sanders et al. 2003). The miss- 
ing galaxies are very similar to those observed, and 
thus this sample is representative of the class of local 
LIRGs as a whole. The sample covers an IR luminos- 
ity range of logii R [ 8 _i oopm] = 11-11.9 [L Q ], and 
includes a variety of morphologies and environments 
(isolated systems, galaxies in groups, and interact- 
ing galaxies), and nuclear activity (Hn, Seyfert, and 
LINER). Details on the observations and near-IR 
properties of the sample are given by Alonso-Herrero 
et al. (2006a, AAH06 hereafter). 

3. OBSERVATIONS 

We obtained T-RcCS imaging observations of 
all the LIRGs in our sample that can be observed 
from Gemini-South. We focus here on the first re- 
sults obtained for four LIRGs (see Alonso-Herrero et 
al. 2006b). Three LIRGs (NGC 5135, IC 4518W, 
and NGC 7130) were observed with the N broad- 
band filter (central wavelength A c = 10.36 /im) and 



NGC 3256 with the Si- 2 narrow-band filter (A c = 
8.74 ism). The T-ReCS observations were effectively 
diffraction limited (FWHM~ 0.30") implying spatial 
resolutions of ~ 50 — 100 pc for our LIRGs. 

In addition we have retrieved archival 
Spitzer/IRAC 8/zm and MIPS 24 /im images 
for our sample, with spatial resolutions of ~ 2" and 
- 5", respectively. The ffST/NICMOS, T-ReCS, 
and Spitzer/IRAC images of three LIRGs are 
presented in Fig. 1. 

4. MID-IR EMISSION ON SCALES OF 
TENS-HUNDREDS PARSECS 
4.1. Morphology 

The morphological resemblance (see Fig. 1) be- 
tween the bright mid-IR regions (traced by T-ReCS) 
and the Hn regions (traced by the NICMOS Paa 
emission) on scales of tens-hundreds of pc in LIRGs 
(see also Soifcr et al. 2001) is remarkable. The 
IRAC 8 /um images, taking into account their lim- 
ited spatial resolution, are similar to the T-ReCs 
images, with the advantage that they are more sen- 
sitive to the diffuse emission not directly associated 
with the bright Hn regions (see §5 for a detailed dis- 
cussion). The near-IR continuum emission is more 
extended than the mid-IR emission, and generally 
does not resemble the morphology of the ionized gas 
(see AAH06). As Ais^m ~ 3.3 x A N , the lack of 
correspondence between the near and mid-IR (and 
Paa) emitting regions is probably not a differential 
extinction effect. Rather, this may reflect age dif- 
ferences since the near-IR continuum is produced 
mostly by stellar populations older than the ioniz- 
ing stellar populations (see e.g., Alonso-Herrero et 
al. 2002), even though highly obscured regions tend 
to be associated with the youngest regions in LIRGs 
(Soifer et al. 2001). 

The nuclei are the brightest mid-IR sources in 
all the galaxies. The Sy2 nuclei of NGC 5135 and 
IC 4518W appear unresolved in the T-ReCs images, 
in agreement with the detection of nuclear point 
sources in the higher spatial resolution (~ 0.15") 
NICMOS images (see AAH06). The nucleus of 
NGC 7130 (classified MS c\ Sy or LINER) appears 
marginally resolved in the T-ReCS images, and the 
NICMOS continuum images reveal the presence of 
at least three sources of similar flux in the central 
~ 190 pc. This suggests that the putative AGN in 
NGC 7130 docs not dominate (see §4.2) the TV-band 
nuclear emission. NGC 3256 is a well-studied merger 
galaxy, with two bright nuclei, which are also the 
brightest mid-IR sources in this galaxy. Although it 
has been argued that this merger may contain ob- 
scured AGN, the extended (FWHM ~ 60 pc) nature 
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Fig. 1. From left to right: HST/NICMOS 1.6 /im or 1.1pm continuum emission, ifST/NICMOS Pa« emission line, 
Gemini/T-ReCS 10.4 /im (TV-band) or 8.7 /im imaging, and Spitzer/IRAC 8 pm imaging of the central regions of three 
LIRGs in our sample. Orientation is north up, east to the left. The two brightest mid-IR and Paa sources of NGC 3256 
are the north and south nuclei of this merger LIRG. For NGC 5135 the bright mid-IR central source is the Seyfert 2 
nucleus, whereas the sources located to the east and south of the nucleus are bright Hn regions (see §5). 



of the near and mid-IR continuum emission of both 
nuclei suggests that such AGN, if present, are not 
dominant in the mid-IR. 

4.2. AGN Contribution to the mid-IR emission 

Approximately 25% of local LIRGs, including 
those in our sample, host spectroscopically confirmed 
AGN (e.g., Veilleux et al. 1995). The Seyfert nucleus 



of IC 4518W appears to be dominate the mid-IR 
emission, whereas for NGC 5135 and NGC 7130 the 
Seyfert nuclei contribute ~ 20 % or less of the N- 
band emission within the central 10". Three more 
galaxies in the LIRG sample of AAH06, not observed 
with T-ReCS, have spectroscopically confirmed Sy 
nuclei. For two of them, the Bl nucleus of the 
IC 694/NGC 3690 (Arp 299) system (Garcfa-Marm 
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Fig. 2. Monochromatic (z/ /„) 8/xm vs. extinction- 
corrected Paa luminosities. The small open and filled 
symbols are the nuclei and Hn regions of LIRGs, re- 
spectively. The large open symbols for the LIRGs in- 
dicate the integrated properties over the NICMOS FOV 
(~ 19" x 19", between 3 and 7kpc for these LIRGs). The 
central 6kpc M51 Hn regions and integrated emission 
from CAL05 are shown as star symbols. The solid line is 
our least-squares fit to the CAL05 M51 H II region data 
extrapolated to the LIRG luminosities. The dashed line 
is the Wu et al. (2005) relation for star-forming galaxies. 

et al. 2006), and NGC 7469, the AGN contribution 
to the mid-IR emission is ~ 30% (Keto et al. 1997; 
Soifer et al. 2003). For the third one, there is no 
mid-IR information. This is consistent with the ma- 
jority of AGN in our local LIRGs not dominating the 
mid-IR emission. 

5. NUCLEAR, HII REGION AND INTEGRATED 
8-10 fjM EMISSION OF LIRGS 

For the nearby galaxy NGC 300 Helou et al. 
(2004) showed that the 8 /im emission tends to be 
more extended than the Ha emission and highlights 
the rims of the Hll regions. They concluded that 
the 8 /im emission is more closely associated with 
photodissociation regions (PDRs) than with ionized 
regions. In Fig. 2 we compare the Paa and 8^m 
emission for resolved Hn regions, nuclei, and inte- 
grated (central 3 — 7 kpc emission of the four LIRGs 
observed with T-ReCs. We also show in this figure 
data for resolved H II knots in M51 from CAL05. The 
LIRG Hn regions, although can be up to 10 times 
more luminous than those in M51, tend to follow the 
extrapolation of the CAL05 relation. Conversely, the 



integrated ~ 3 — 7 kpc emission of LIRGs deviates 
significantly from the relation found for individual 
Hn regions, but only slightly from the fit of Wu et 
al. (2005) for star- forming galaxies (Fig. 3). This 
could arise from two causes. First, Wu et al. (2005) 
may have underestimated the reddening (obtained 
from the Balmer decrements) in these dusty star- 
forming galaxies. Second, the 8 /xm emission may be 
dominated by diffuse emission not associated directly 
with the H II regions (see below) . 

The LIRG nuclei (except the south nucleus of 
NGC 3256) show elevated mid-IR/Paa ratios when 
compared to Hll regions. Such spatial differences 
in the mid-IR emission of nuclear and circumnuclear 
regions are observed for our LIRGs from high-spatial 
resolution T-ReCS spectroscopy (sec Diaz-Santos et 
al. 2006). One possibility is that an insufficient ex- 
tinction correction (as nuclear extinctions in LIRGs 
tend to be higher than in extranuclear regions) will 
produce a differential effect, making the nuclei ap- 
pear more IR-luminous. The AGN present in three 
of our LIRGs may also play a role, as their continua 
are produced by dust heated to higher temperatures 
than Hn regions, and do not present the strong PAH 
emission characteristic of Hn+PDRs (e.g., Roche et 
al. 2006; Diaz-Santos et al. 2006). 

Since a variety of Hn regions in M51 and LIRGs 
have a similar 8 /jm vs. Paa relation we conclude 
that their 8 (im emission in the Spitzer/IRAC band is 
well characterized by a thermal continuum plus PAH 
features with no strong variations over the range of 
conditions probed here (e.g., metalicity near or over 
solar, see AAH06 and CAL05). 

However, the integrated central 3 — 7 kpc 8 /im 
vs. Paa emission differs significantly in all these 
environments from the relation found for the indi- 
vidual Hn regions. This may be explained by the 
presence, in addition to the bright and compact H II 
regions, of a more diffuse and extended 8 /im compo- 
nent (see Fig. 1, and Helou et al. 2004; Diaz-Santos 
et al. 2006). This extra emission at 8/iin would 
be produced not by local, strong ionizing sources, 
but by the diffuse radiation field that permeates the 
ISM. As such, the spectrum in these regions would 
be characterized by a weak continuum and strong 
PAH features with a large equivalent width. This 
has already been confirmed for one of the LIRGs in 
our sample as discussed in detail by Diaz-Santos et 
al. (2006). Hence when compared to individual Hn 
regions, an excess of 8 fim/Paa emission can be ex- 
pected for the integrated properties over a few kpc. 
This is further supported by the fact that the cen- 
tral 3 — 7 kpc emission of LIRGs falls only slightly 
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Fig. 3. Monochromatic 24 (jm luminosities vs. extinction-corrected Paa luminosities (left) for LIRGs, normal galaxies 
(Boker et al. 1999), ULIRGs (Murphy et al. 2001) ULIRGs, and Hn regions in M51 and M81 (see text and AAH06). 
The solid line is the fit to the M51 Hn and the normal galaxies and LIRGs, and dashed line the CAL05 fit to the M51 
Hn knots. The 24^tm/Paa ratio vs. the extinction-corrected Paa luminosity (right). Symbols are as in left panel. 



below the Ha vs. 8 relation found for the in- 
tegrated properties of the galaxies studied by Wu 
et al. (2005), as indicated in Fig. 2. We conclude 
that calibrations of the SFR in terms of the 8 /im 
emission for distant galaxies should be based on the 
integrated mid-IR emission of galaxies, not that of 
the H II regions alone. 

6. THE 24 EMISSION AS A SFR TRACER 
FOR DUSTY STAR-FORMING GALAXIES 

Hclou et al. (2004) showed that in contrast 
with the 8 /um emission, the 24 /im emission of 
NGC 300 appears to be more centrally peaked in 
star-forming regions. In Fig. 3 (left) we com- 
pare the monochromatic 24 /im luminosities and the 
extinction-corrected Paa luminosities (see AAH06) 
for all the LIRGs in our sample. We also include 
a small comparison sample of normal galaxies and 
ULIRGs, and data for resolved star-forming regions 
in M51 from CAL05, and in M81 from PGPG06. 

The integrated emission of LIRGs, ULIRGs, and 
normal galaxies seems to continue the relation be- 
tween the extinction-corrected Paa luminosity and 
the 24 /xm luminosity observed by CAL05 for the 
M51 central Hn knots, unlike the case of the 8/zm 
emission. The fit to the LIRGs, ULIRGs, normal 
galaxies, and M51 Hn knots provides a slope of 
1.148 ±0.013 (Fig. 3; Wu et al. 2005; AAH06), only 
slightly larger than the slope of the CAL05 fit to 



the M51 Hn knots. The M81 Hn regions, on the 
other hand, appear to be offset with respect to M51 
Hn regions, normal galaxies, and LIRGs. This be- 
havior could arise from a lower UV absorption ef- 
ficiency in the relatively low luminosity and lightly 
obscured Hn regions in M81, as well as to uncertain 
extinction corrections (see discussion by PGPG06). 
There are galaxy-to-galaxy variations in the mid-IR 
to Paa ratios (Fig. 3, right), as also noted by CAL05. 
The variation of the £(24/zm)/L(Paa) corr ratio with 
the Paa luminosity (i.e., SFR) between the M51 Hn 
knots and our LIRGs and ULIRGs is similar to the 
^24^im/SFR ratios found by CAL05. In fact, the gap 
between the M51 Hn knots and the LIRGs/ULIRGs 
seen in our figure is populated by their UV-selected 
starbursts. 

There are two possibilities which might explain 
the deviation from strict proportionally of the 24 fj.m 
vs. Paa relation. First the nuclear extinctions of 
LIRGs may be underestimated (the values averaged 
over the Paa emitting regions are Ay ^2 — 6 mag, 
see AAH06), whereas some LIRGs and ULIRGs are 
known to contain highly obscured regions, usually as- 
sociated with the nuclei of the galaxies (e.g., Gcnzcl 
et 1998; Alonso-Hcrrcro et al. 2000). From the ob- 
served ratios of mid-IR emission lines of star-forming 
galaxies there is evidence that some of the youngest 
stars in LIRGs may still embedded in high density ul- 
tracompact H n regions and hidden from us by large 



6 



ALONSO-HERRERO 



amounts of extinction (Rigby & Rieke 2004). In such 
regions the dust competes increasingly effectively for 
ionizing photons and UV continuum photons, so that 
an increasing fraction of the luminosity is expected 
to emerge in the IR. The second possibility would be 
if the LIRG emission is dominated by AGN, but as 
discussed in §4.2 AGN do not appear to dominate 
the mid-IR emission in our sample of LIRGs. 

6.1. The empirical calibration of SFR vs. 24 fim 

In deriving his relationship between the SFR and 
the IR luminosity Kennicutt (1998) assumed that the 
great majority of the luminosity from young stars 
would be absorbed by dust and reradiated in the far- 
IR. Although this assumption is likely to be correct, 
observationally there may be other contributions to 
the total IR luminosity from older populations of 
stars or other luminosity sources. Therefore, a true 
total IR luminosity measurement may overestimate 
the recent star formation in a galaxy, although es- 
timates based on IRAS measurements alone are less 
subject to these issues because they only poorly sam- 
ple the output of the cold dust (because the longest 
band is at 100/im). Second, Kennicutt assumed a di- 
rect proportionality between the Ha luminosity and 
the total IR. We find empirically that the increasing 
absorption efficiency in more luminous and obscured 
galaxies leads to a deviation from strict proportion- 
ality, causing the IR output to rise with increasing 
star-forming luminosity. 

Assuming Case B recombination, and using the 
relation between SFR and Ha quoted by Kennicutt 
(1998) and the fit between the 24/itm and Paa lu- 
minosities (AAH06) we derive the following relation 
between the SFR rate and the 24 fan luminosity for 
luminous, dusty galaxies: 

SFR = 8.45 x 10" 38 (£(24 fim)) 871 (1) 

where the SFR is in Mgyr -1 and the 24 /Urn lumi- 
nosity is in ergs -1 . This relation is analogous to 
the widely used relation between SFR and IR lumi- 
nosity (Kennicutt 1998) but it is not affected by the 
uncertain contribution to the total IR luminosity of 
a galaxy of dust heating from old stars. 

Summarizing, we have used the sensitivity of 
Spitzer and the high spatial resolution of T-ReCS 
on the Gemini-South telescope to study the mid-IR 
properties of LIRGs. The T-ReCS 8 — 10 /im imaging 
observations of LIRGs have allowed us to resolve the 
nuclear emission (star formation and/or AGN) from 
that of circumnuclear Hn regions. The comparison 
of the 8 /rni/Paa ratios of the integrated and resolved 



H II regions of LIRGs reveals an 8 [im diffuse compo- 
nent, not directly related to the ionizing stars, that 
can be as luminous as that from the resolved Hn 
regions. Therefore the calibration of the SFR for 
distant galaxies should be based on the integrated 
mid-IR emission of nearby galaxies, not that of H II 
regions alone. We provide a calibration of the SFR 
in terms of the integrated 24 ftm luminosity. Simi- 
lar mid-IR studies with the GTC/CanariCam system 
will be possible in the near future. 
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